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(54) Method of calibrating a logging tool 

(57) A logging tool that is exposed to an external en- 
ergy spectrum is calibrated by (a) providing a stabiliza- 
tion source that emits a stabilization signal having a 
known energy value; (b) providing a receiver that re- 
ceives a total spectrum that comprises the sum of the 
external energy spectrum and the received stabilization 
signal; (c) adjusting the gain of the tool so that the re- 



ceived stabilization signal has said known energy value, 
by the steps of: (ci) ensuring that the received stabiliza- 
tion signal is within a preferred range of said known en- 
ergy value; (cii) calculating a value D n that is proportion- 
al to a derivative of the spectrum at the stabilization peak 
energy, where n is the order of the derivative and is odd; 
(ciii) adjusting the gain of the tool on the basis of the 
value of D- 
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Description 

[0001] The present invention relates generally to a method of calibrating a logging tool, particularly but not exclusively 
a gamma density logging tool. 

5 [0002] Gamma density logging is well known in the art of petroleum drilling. In gamma density logging, a source of 
gamma rays is lowered into the hole, along with at least one, and more typically at least two gamma ray detectors. 
Each detector typically comprises a large thallium-doped sodium iodide NaI(T1) crystal coupled to a low-noise photo- 
multiplier. The incidence of gamma rays on the detectors, along with known information about the respective distances 
of each detector from the GR source, gives information about the effect of the surrounding formation on the gamma 

10 rays as they are transmitted from the source to the detectors, and thus gives information about the formation itself. 
[0003] In nuclear radiation measurement techniques and also in photometry use is made of detectors such as scin- 
tillation counters, proportionally counter tubes, semiconductors, or photomultipliers having internal amplification, which 
supply pulses whose pulse amplitude corresponds to the energy absorbed by the detector. Gamma density logging 
tools typically use a gamma ray source in conjunction with one or more scintillation counters. Gamma rays from the 

is source are emitted into the formation surrounding the tool. The formation scatters the gamma rays, causing some of 
them to return to the tool and impact the scintillation counters. The counters convert energy from these incident gamma 
rays into light, which is in turn detected and amplified by a photomultiplier, which emits a signal corresponding to the 
energy of the incident gamma rays. 

[0004] A particular difficulty arises as a result of drift in such radiation measurement means, especially because of 

20 variations in internal amplification. These variations are particularly severe as the temperature of the tool rises. These 
variations gives rise to variations in pulse heights and consequently to a displacement of the energy lines or energy 
peaks with respect to the adjusted response thresholds of pulse registration, and hence to inaccurate measurements 
of intensity. The variations can be kept within acceptable limits by regular checking with calibrating radiation sources 
and adjustment of the tool. However, periodic calibration of the energy spectrum with external sources is not possible 

25 when the tool is in use (logging). For this reason, and because the gain changes due to elevated temperatures in a 
well can be quite large, automatic control methods and devices for stabilization of drift are preferred. 
[0005] In order to avoid the disadvantages associated with drift, it is desirable to use a peak of the pulse amplitude 
spectrum as the controlling variable in automatic stabilization techniques. For this to be possible, the pulse amplitude 
spectrum of the radiation must have a peak that is discemable and suitable for control purposes. 

30 [0006] Typical gamma-ray logging tools use a cesium-1 37 gamma source and two density detectors spaced apart 
from the source. Cesium-1 37 emits gamma rays with an energy of 662 keV. The standard means for stabilizing such 
a detector in a gamma-gamma density tool is to place a small cesium-1 37 source on or near the surface of each 
detector. A representative spectrum of the stabilization source energies detected at the far detector is shown in Figure 
1 . It can be seen that there is a clear peak at 662 keV, but there are also a large number of counts at lower energy 

35 levels, which correspond to smaller amounts of energy that were absorbed in the scintillation crystal. 

[0007] In operation, the observed spectrum for a tool of this sort comprises the sum of the logging and stabilization 
spectra. Fortunately, the energy loss in the formation is great enough that the energy of the gamma rays detected after 
passage through the formation is much lower than their original energy, as shown in Figure 2. In contrast, the energy 
of the gamma rays originating at the stabilization source does not pass through the formation and so is virtually un- 

^0 changed when it is detected. The observed spectrum thus includes a stabilization peak at 662 keV superimposed on 
a spectrum having a much larger peak at lower energies, as in the spectrum of Figure 3, in which the portion of the 
spectrum between 40 and 400 keV is off the scale. 

[0008] It is presently known to use a two step process for distinguishing the stabilization peak from the rest of the 
spectrum. The two steps entail (1 ) adjusting the gain of the system so as to place the stabilization peak at approximately 
45 662 keV and (2) adjusting the gain until the count rate in the energy range from 617 to 662 keV is the same as the 
count rate in the energy range from 662 to 708 keV Step (1 ) is carried out by positioning the peak such that there are 
enough counts above 617 keV to ensure that the stabilization peak is above that point and so few counts above 708 
to ensure that the peak is below that point. 

[0009] Successful use of this technique depends on three factors. First, there must be relatively few counts above 
so the stabilization peak. Second there must be a well-defined valley between the logging spectrum and the stabilization 
peak, and third, the stabilization peak must be easily distinguished in the spectrum. In theory, these three conditions 
can always be met by increasing the strength of the stabilization source. Increasing the strength of the stabilization 
source also increases the fraction of the spectrum that corresponds to the stabilization source, however. Since the 
stabilization spectrum extends into the energy region that is used to compute formation density, the portion of the 
55 spectrum that is attributable to the stabilization source must be subtracted from the observed spectrum to allow cal- 
culation of the formation density. Thus, as the strength of the stabilization source increases, the statistical uncertainty 
of the desired calculations also increases. This is particularly disadvantageous in situations where the overall spectrum 
has a low count rate, such as in logging while drilling (LWD) applications and in cased-hole density tools. 
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[001 0] Cased-hole density measurements are particularly difficult, since the tool must operate in both very high- and 
very low-count-rate situations. This is because gamma rays are attenuated exponentially with the density of the rock. 
High-density formations, such as carbonates, attenuate gamma rays quickly. However, very low density formations, 
such as coal, attenuate the gamma rays much less. Thus, if a detector far from the source is large enough to detect a 

5 reasonable number of counts per second in a dense rock, it will detect many times more counts in a low-density rock. 
If the tool is used in open holes without casing, the fractional difference between high-density formations behind casing 
and low-density formations in open holes is even greater, since the casing also causes significant attenuation. Even 
if the tool is not used in open holes, it is desirable to calibrate it in existing calibration fixtures, which have no casing. 
Hence, the tool must be able to operate under a wide range of conditions. 

10 [0011] Furthermore, if the tool uses a cobalt-60 logging source, which emits gamma rays of 1.17 and 1.33 MeV, the 
high -count-rate situations have a significant number of counts at energies in the region above 662 keV. Figure 4, which 
is an example of an observed spectrum for the third detector of a cased-hole density tool placed in a magnesium block, 
illustrates the nature of this problem. As can be seen, there are a significant number of counts above 662 keV and 
there is no well-defined valley between the stabilization peak and the logging spectrum. In situations like this, the 

is conventional two-step technique described above for adjusting gain will not work. In addition, because the stabilization 
peak is obscured by a large, nonlinear background, the technique of balancing the count rates immediately above and 
below 662 keV will result in a miscalibrated instrument. 

[0012] Hence, it is desired to provide a stabilization technique that is effective in situations where there are a signif- 
icant number of counts above 662 keV, where there is no well-defined valley between the stabilization peak and the 
20 logging spectrum, and/or where the stabilization peak is obscured by a large, nonlinear background. The stabilization 
technique should require sufficiently few calculations as to be operable in the logging context and should be consistently 
accurate in a variety of logging situations. 

[001 3] We have now devised a method of discerning and calibrating the stabilization peak in a gamma density logging 
tool, which method provides the advantage of being equally effective in both low and high count rate situations. Another 
25 advantage of the present method is that it accurately locates the stabilization peak even in situations where there are 
significant amounts of energy above the stabilization peak, where no well-defined valley sets off the stabilization peak, 
and where a large non-linear background obscures the stabilization peak. 

[0014] In one aspect, the invention provides a method of calibrating a logging tool that is exposed to an external 
energy spectrum, comprising the steps of: (a) providing a stabilization source that emits a stabilization signal having 

30 a known energy value; (b) providing a receiver that receives a total spectrum that comprises the sum of the external 
energy spectrum and the received stabilization signal; (c) adjusting the gain of the tool so that the received stabilization 
signal has said known energy value, by the steps of: (ci) ensuring that the received stabilization signal is within a 
preferred range of said known energy value; (cii) calcu lating a value D n that is proportional to a derivative of the spectrum 
at the stabilization peak energy, where n is the order of the derivative and is odd; (ciii) adjusting the gain of the tool on 

35 the basis of the value of D,,. 

[0015] In another aspect, the invention provides a method of calibrating a logging tool that is exposed to an external 
energy spectrum, comprising the steps of: (a) providing a stabilization source that emits a stabilization signal having 
a known energy value; (b) providing a receiver that receives a total spectrum that comprises the sum of the external 
energy spectrum and the received stabilization signal; (c) adjusting the gain of the tool so that the received stabilization 

40 signal has said known energy value, by the steps of: (ci) calculating the centroid of the area under the total spectrum; 
(cii) adjusting the gain of the tool on the basis of the calculated centroid; (ciii) dividing a portion of the spectrum that 
includes the known energy value of said stabilization source into a number of windows of equal width, with an equal 
number of windows on either side of the known energy value; (civ) using the count rates in the windows established 
in step (ciii) to calculate a value D„ that is proportion to the derivative of the spectrum at the stabilization peak 

45 energy, n preferably being 3; (cv) adjusting the gain of the tool on the basis of the value of In a further aspect, the 
invention provides a method of calibrating a logging tool that is exposed to an external energy spectrum, comprising 
the steps of: (a) providing a stabilization source that emits a stabilization signal having a known energy value; (b) 
providing a receiver that receives a total spectrum that comprises the sum of the external energy spectrum and the 
received stabilization signal; (c) adjusting the gain of the tool so that the received stabilization signal has said known 

50 energy value, by the steps of: (ci) ensuring that the received stabilization signal is within a preferred range of said 
known value; (cii) dividing a portion of the spectrum that includes the known energy level of the stabilization peak into 
a number of windows of equal width, with an equal number of windows on either side of the stabilization peak energy; 
(ciii) using the count rates in the windows established in step (cii) to calculate a value D n that is proportional to the n m 
derivative of the spectrum at the stabilization peak energy; (civ) adjusting the gain of the tool on the basis of the value 

55 of and (cv) repeating steps (cii) - (civ) until D n equals substantially zero. 

[001 6] In order that the invention may be more fully understood, reference will hereafter be made to the accompanying 
drawings wherein: 
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Figure 1 is a typical spectrum of energy from a stabilization source placed on the far detector of a wireline density 
tool; 

Figure 2 is a logging spectrum from the far detector of a wireline density tool that is placed in a magnesium block; 
Figure 3 is a combined logging and stabilization spectrum from the far detector of a density tool in a magnesium 
s block; 

Figure 4 is a an example of an observed (combined) spectrum for the third detector of a cased-hole density tool 
placed in a magnesium block; and 

Figure 5 is an annotated combined spectrum illustrating the technique of the present invention. 

10 [001 7] The present invention comprises a technique for processing logging data from a stabilized source. Thus, the 
input data may resemble the spectra shown in Figures 3 and 4, discussed above. Two objectives of the technique are 
to discern the stabilization peak and to accurately position it at 662 keV The present technique uses two steps, both 
of which are different from the steps of previously known stabilization techniques. 

[0018] The first step of the present technique entails positioning the peak at approximately the correct position by 
is computing the centroid of the area under the entire spectrum curve. This is effective because the shape of the observed 
logging spectrum is relatively insensitive to the density and lithology of the formation, so that the spectra of various 
formations should have centroids at approximately the same energy. 

[0019] The centroid of the spectrum is calculated by first dividing the spectrum into a plurality of contiguous energy 
windows that together span most of the spectrum. While it is preferable to use a large number of windows that span 

20 the entire spectrum, this must be balanced against the advantages of simplifying the necessary computation. In view 
of these competing considerations, it has been found that a preferred number of windows is between 8 and 16. The 
highest-energy window preferably encompasses all energies above its threshold. The lower end of the lowest-energy 
window must be set above the noise threshold for that detector and its expected operating range. This is typically about 
40 keV. An example of a suitable division is illustrated in Figure 5. 

25 [0020] Once the windows have been set, the centroid is calculated. According to a preferred embodiment, the centroid 
is calculated according to the following equation: 



so Y*y* x * 

centroid = 

i 

35 

where x- is the energy of window / and y- is the count rate in that window. The calculated centroid is compared to a 
predicted nominal value, and the difference is used to adjust the gain of the instrument. The nominal position for the 
centroid depends on the detector geometry and tool design. It is determined by acquiring data in a wide range of 
formations, adjusting the gain of the measured spectrum so that the location of the stabilization peak is at the desired 
40 value, and then computing the centroid. The range of centroids determined for the various formations is then used to 
generate a nominal centroid position for a given formation. 

[0021] According to one embodiment, a preferred first step is to determine whether the stabilization peak is close 
enough to the desired (nominal) position to allow the first (coarse) tuning step to be omitted, in which case the third 
derivative can be used to fine-tune the stabilization. This determination can be made by examining the count rates in 

45 the windows Stabl and above (Figure 5). For instance, if the Gamma window has more than some minimum fraction 
of the total counts, if Stab2 + Stab3 have many fewer counts than are expected in the stabilization peak, or if Stab4 
counts are greater than those in Stab3, then the peak is clearly not close enough to its target value. On the other hand, 
if Stab 1 is less than Stab2 or Stab2 is less than 0.92 x Stab3, then the peak is close enough to be fine tuned, as 
described below. If the peak is not close enough, the centroid is first adjusted in a coarse tuning step as follows. 

so [0022] The centroid is adjusted linearly in proportion to the ratio of the difference between the calculated and nominal 
values of the centroid to the calculated centroid value. For example, if the calculated centroid occurs at 200 keV and 
the nominal value is 250 keV, it is necessary to increase the gain by 25 percent [(250-200)/200]. Because the centroid 
calculation is not precise, this process of calculating the centroid and adjusting the gain is preferably iterated until the 
difference between the calculated centroid and the nominal value is small, such as on the order of four percent (4%) 

55 of the nominal value. Alternatively, a predefined narrow target range that includes the nominal value can be used, with 
the gain being adjusted until the calculated centroid falls within this range. By way of example only, the target range 
can be defined as the nominal value plus or minus 5 keV. A narrower target range, such as the nominal value plus or 
minus 2 keV could alternatively be used. 
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[0023] After this step is completed, it may be necessary to make additional minor adjustments to account for small 
variations in the shape of the spectrum prior to completion of the coarse adjustment step. One of two types of adjustment 
is preferably used. One type of adjustment applies in the case where the measured spectrum includes only the stabi- 
lization source and no logging source is present. This condition occurs when the background spectrum is being meas- 
5 ured during the calibration procedure. If the total count rate is so low that no logging source is likely to be present, the 
spectrum is shifted an additional amount, based on a knowledge of where the centroid for the stabilization peak should 
be. 

[0024] A second type of adjustment may be required to account for spectral shifts due to the lithology of the formation. 
This adjustment can be based on the ratio of the count rates in a pre-selected low-energy window to a pre-selected 
10 high-energy window. For example, the ratio of the total count rate in the window labeled Den Hi in Figure 5 to the total 
count rate in the window labeled DenLo can be used. If the ratio is lower than a predetermined nominal value, the gain 
should be increased to compensate for the shape of the spectrum and, conversely, if the ratio is greater than a prede- 
termined nominal value, the gain should be decreased. 

[0025] Since the coarse tuning step uses the centroid of the peak to approximate the peak position, it does not require 
is either small count rates above 708 keV or a well-defined valley between the stabilization peak and the bulk of the 
logging spectrum. 

[0026] Once the gain has been adjusted coarsely in this manner, the windows labeled Stabl to Stab4 in Figure 5 
are used to fine tune the stabilization. The windows Stabl through Stab4 are preferably selected so that they are 
contiguous, of equal width, and positioned so that Stabl and Stab 2 are on one side of 662 keV and Stab 3 and Stab 

20 4 are on the opposite side of 662 keV, as shown in Figure 5. A preferred width for each of windows Stabl through 
Stab4 is equal to from about 0.5 to about 1 .5 times the width of the stabilization peak at one-half of its maximum value 
("full width at half max," hereinafter *W). The width of the windows is preferably about 1 .25 times the full width at half 
max (1 .25 W). This quantity can be generated from measurements of the stabilization peak alone, in the absence of a 
logging signal, such as is shown in Figure 1. Alternatively, in many instances, the width of windows Stab 1 through 

25 Stab4 can be set quantitatively, such as for example, in the range of from about 30 keV to about 80 ke V, with a preferred 
width being approximately 45 ± 5 keV. 

[0027] According to a preferred embodiment, the count rates from these windows are preferably combined to yield 
a quantity D 3 that is proportional to the third derivative. While any analysis of the count rates in the Stab windows that 
yields a value for D 3 that is proportional to the third derivative can be used, a preferred equation for calculating D 3 is 
30 as follows: 

D 3 - Stabl - 3Stab2 + 3Stab3 - Stab4, 

35 Since the stabilization peak is symmetric, its third derivative at 662 keV is zero. Thus, the gain is adjusted until D 3 
equals zero. Because the third derivative of a quadratic is also zero, this technique will place the center of the stabili- 
zation peak at 662 keV, even if it is superimposed on a relatively large quadratic background. Hence, even relatively 
small stabilization peaks are discernable and useable, with the result that tools using the present technique can operate 
with much smaller stabilization sources than were previously useful. 

40 [0028] Alternatively, the mathematical techniques for fine tuning the tool can be varied without departing from the 
concepts of the present invention. For example, other derivatives can be used. It is preferred that odd-ordered deriv- 
atives be used, so the derivative will equal zero at the stabilization peak. Also, an even number of windows is required 
to compute an odd derivative. 

[0029] For example, using the first derivative requires only two windows, Stabl and Stab2, and gives a value for D 1 
45 according to the equation: 

D 1 = Stabl - Stab2 

so Using this analysis, D1 will be zero at the peak only if there is no background noise. Because there is typically some 
level of background noise, use of the first derivative is not preferred. 

[0030] Higher-order derivatives can also be used, but they require more windows and the statistical uncertainty of 
the measurement increases. The advantage of higher orders is that the higher the order, the more rapidly the back- 
ground can vary without affecting the measurement of the derivative. Examples of higher-order equations are: 



55 



D 6 = Stabl - 5Stab2 + 10Stab3 - 10Stab4 + 5Stab5 - Stab6 and 
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D 7 = Stabl - 7Stab2 + 21 Stab3 - 35Stab4 + 35Stab5 - 21 Stab6 + 7Stab7 - Stab8 

It will be understood by those skilled in the art that these equations are merely representative variations and embodi- 
5 ments of the mathematical technique and are not intended to limit the invention to any particular embodiment. 

[0031] As discussed above, the present method is effective for discerning the stabilization peak even when there 
are high count rates above the stabilization peak, or where there is no well-defined valley between the logging spectrum 
and the stabilization peak. These conditions can be determined objectively, as for example when the number or fraction 
of counts above a certain threshold exceeds a preset limit, or, preferably, subjectively on the basis of the spectrum 
io itself. For example, for purposes of the present discussion, the condition of high count rates above the stabilization 
peak can be defined as occurring when the amplitude of the spectrum at a point 2 W above the stabilization peak is 
greater than 20% of the maximum amplitude of the peak due only to the stabilization source. Likewise, for purposes 
of the present discussion, lack of a well-defined valley occurs when no portion of the spectrum between the stabilization 
peak and the logging peak is less than 40% of the maximum amplitude of the peak due only to the stabilization source. 
is it will be understood that these conditions are merely descriptive of conditions under which the present technique will 
work; it is not necessary or preferred to determine their presence in actual operations. 

[0032] It will be understood by those skilled in the art that it is preferred that the processes described be carried out 
iteratively and continuously during logging, such as, for example, every half-second. The frequency of the calculation 
can be varied, but is preferably at least once per minute. 
20 [0033] While preferred embodiments of the present method have been disclosed herein, ft will be understood that 
modifications thereof can be made by one skilled in the art without departing from the invention as defined in the claims. 



Claims 

25 

1 . A method of calibrating a logging tool that is exposed to an external energy spectrum, comprising the steps of: (a) 
providing a stabilization source that emits a stabilization signal having a known energy value; (b) providing a re- 
ceiver that receives a total spectrum that comprises the sum of the external energy spectrum and the received 
stabilization signal; (c) adjusting the gain of the tool so that the received stabilization signal has said known energy 
30 value, by the steps of: (ci) ensuring that the received stabilization signal is within a preferred range of said known 

energy value; (cii) calculating a value D n that is proportional to a derivative of the spectrum at the stabilization 
peak energy, where n is the order of the derivative and is odd; (ciii) adjusting the gain of the tool on the basis of 
the value of D n . 

35 2. A method according to claim 1 , wherein step (ci) is carried out by (ci)[1 ] dividing said total spectrum into a plurality 
of adjacent windows, and (ci)[2] evaluating the count rates in said windows. 

3. A method according to claim 2, wherein step (ci)[2] comprises either comparing the count rate in the high end of 
the spectrum to the count rate of the total spectrum; or comparing the count rate within a predetermined range of 

40 a predicted stabilization peak to a predetermined expected value for said count rate; or comparing the count rate 

in a window Stab3 that is adjacent to and above a predicted stabilization energy to the count rate in a window 
Stab4 that is adjacent to and above Stab3; or comparing the count rate in a window Stab2 that is adjacent to and 
below a predicted stabilization energy to the count rate in a window Stabl that is adjacent to and below Stab2; or 
comparing the count rate in a window Stab2 that is adjacent to and below a predicted stabilization energy to the 

45 count rate in a window Stab3 that is adjacent to and above said predicted stabilization energy. 

4. A method according to claim 1 , wherein step (ci) is carried out by: (ci)[a] calculating the centroid of the area under 
the total spectrum; and (ci)[b] adjusting the gain of the tool on the basis of the calculated centroid. 

50 5. The method according to claim 4, wherein the centroid of the area under the total spectrum is calculated by defining 
a plurality of contiguous energy windows that together span a majority of the spectrum and calculating the centroid 
according to the following equation: 



55 
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centroid = -i= 

2> 



where x- is the energy of window / and y/ t is the count rate in that window. 

10 6. A method according to claim 1 , wherein step (cii) is carried out by dividing a portion of the spectrum that includes 
the known energy level of the stabilization peak into a number, preferably (n+1), of windows of equal width, with 
an equal number of windows on either side of the stabilization peak energy and calculating a value for D n using 
the count rates in each of the windows. 



is 7. A method according to claim 1, wherein n = 3 and a value for D 3 is calculated according to the equation: D 3 = 
Stabl - 3Stab2 + 3Stab3 - Stab4, where Stabl through Stab4 are consecutive, contiguous windows of equal width 
defined such that the known stabilization peak energy lies between Stab2 and Stab3; or n = 5 and a value for D 5 
is calculated according to the equation: D 5 = Stabl - 5Stab2 + 10Stab3 - 10Stab4 + 5Stab5 - Stab6, where Stabl 
through Stab6 are consecutive, contiguous windows of equal width defined such that the known stabilization peak 

20 energy lies between Stab3 and Stab4; or n = 7 and a value for D; is calculated a according to the equation: D? = 

Stabl - 7Stab2 + 21 Stab3-35Stab4 + 35Stab5 - 21Stab6 + 7Stab7 - Stab8, where Stabl through Stab8 are con- 
secutive, contiguous windows of equal width defined such that the known stabilization peak energy lies between 
Stab4 and StabS. 



2S 8. A method of calibrating a logging tool that is exposed to an external energy spectrum, comprising the steps of: (a) 
providing a stabilization source that emits a stabilization signal having a known energy value; (b) providing a re- 
ceiver that receives a total spectrum that comprises the sum of the external energy spectrum and the received 
stabilization signal; (c) adjusting the gain of the tool so that the received stabilization signal has said known energy 
value, by the steps of: (ci) calculating the centroid of the area under the total spectrum; (cii) adjusting the gain of 

30 the tool on the basis of the calculated centroid; (ciii) dividing a portion of the spectrum that includes the known 

energy value of said stabilization source into a number of windows of equal width, with an equal number of windows 
on either side of the known energy value; (civ) using the count rates in the windows established in step (ciii) to 
calculate a value D„ that is proportion to the rt* 1 derivative of the spectrum at the stabilization peak energy, n 
preferably being 3; (cv) adjusting the gain of the tool on the basis of the value of 0^ 

35 

9. A method according to claim 8, wherein the centroid of the area under the total spectrum is calculated by defining 
a plurality of contiguous energy windows that together span a majority of the spectrum and calculating the centroid 
according to the following equation: 

40 

centroid = 

45 t 



where x b is the energy of window / and y,- is the count rate in that window 

so 1 0. A method according to claim 8, wherein the gain is adjusted in step (cii) by an amount proportional to the difference 
between the calculated centroid and a predicted nominal value for said centroid. 

11. A method according to claim 8, wherein steps (ci) and (cii) are repeated until the centroid falls within a predeter- 
mined range that includes a predicted nominal value for said centroid. 

55 



7 



EP 0 942 296 A2 



centroid = 



2>, 



12. A method according to claim 8, further including the step of compensating for the lithology of the formation by 
calculating the ratio of the count rates in a pre-selected low-energy window to the count rate in a pre-selected 
high-energy window, and adjusting the gain on the basis of said ratio. 

13. A method according to claim 8, wherein n = 3 and D 3 is calculated according to the equation: D 3 = Stabl - 3Stab2 
+ 3Stab3 - Stab4, where Stabl through Stab4 are consecutive, contiguous windows of equal width defined such 
that the known stabilization peak energy lies between Stab2 and Stab3; or n = 5 and a value for D 5 is calculated 
according to the equation: D 5 = Stabl - 5Stab2 + 10Stab3 - 10Stab4 + 5Stab5 - Stab6, where Stabl through Stab6 
are consecutive, contiguous windows of equal width defined such that the known stabilization peak energy lies 
between Stab3 and Stab4; or n = 7 and a value for Dj is calculated a according to the equation: Dy = Stabl - 
7Stab2 + 21 Stab3 - 35Stab4 + 35Stab5-21 Stab6 + 7Stab7 - StabS, where Stabl through StabS are consecutive, 
contiguous windows of equal width defined such that the known stabilization peak energy lies between Stab4 and 
StabS. 

14. A method according to claim 8, wherein steps (ciii) - (cv) are repeated until D„ equals substantially zero. 

15. A method according to claim 8, wherein the total spectrum includes high count rates above the stabilization peak; 
and/or the stabilization peak is superimposed on a large quadratic background; and/or there is no well-defined 
valley between the logging spectrum and the stabilization peak. 

16. A method of calibrating a logging tool that is exposed to an external energy spectrum, comprising the steps of: (a) 
providing a stabilization source that emits a stabilization signal having a known energy value; (b) providing a re- 
ceiver that receives a total spectrum that comprises the sum of the external energy spectrum and the received 
stabilization signal; (c) adjusting the gain of the tool so that the received stabilization signal has said known energy 
value, by the steps of: (ci) ensuring that the received stabilization signal is within a preferred range of said known 
value; (cii) dividing a portion of the spectrum that includes the known energy level of the stabilization peak into a 
number of windows of equal width, with an equal number of windows on either side of the stabilization peak energy; 
(ciii) using the count rates in the windows established in step (cii) to calculate a value D n that is proportional to the 
^derivative of the spectrum at the stabilization peak energy; (civ) adjusting the gain of the tool on the basis of 
the value of D„; and (cv) repeating steps (cii) - (civ) until D n equals substantially zero. 

17. A method according to claim 16, wherein step (ci) is accomplished by calculating the centroid of the area under 
the total spectrum and adjusting the gain of the tool on the basis of said calculated centroid, said centroid being 
calculated by defining a plurality of contiguous energy windows that together span a majority of the spectrum and 
calculating the centroid according to the following equation: 



where Xj is the energy of window / and y) is the count rate in that window. 

18. A method according to claim B, further including the step of compensating for the lithology of the formation by 
calculating the ratio of the count rates in a pre-selected low-energy window to the count rate in a pre-selected 
high-energy window, and adjusting the gain on the basis of said ratio. 



centroid = -~ 
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Figure I . Spectrum from a stabilization source placed on the far detec 
tor of a wireline density tool. 
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Figure 2. Logging spectrum from the far detector of a wireline density 
tool that is placed in a magnesium block. The effects of the stabilization 
source have been removed. 
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Figure 3. Combined logging spectrum and stabilization spectrum from 
the far detector of a density tool in a magnesium block. The peak at 662 
is from the stabilization source. The spectrum from 40 to 400 keV is off 
scale. 
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Figure 4 Combined Logging and Stabilization spectrum for the third 
detector of a cased-hole density tool placed in a magnesium block. 
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(54) Method of calibrating a logging tool 

(57) A logging tool that is exposed to an external en- 
ergy spectrum is calibrated by (a) providing a stabiliza- 
tion source that emits a stabilization signal having a 
known energy value; (b) providing a receiver that re- 
ceives a total spectrum that comprises the sum of the 
external energy spectrum and the received stabilization 
signal; (c) adjusting the gain of the tool so that the re- 
ceived stabilization signal has said known energy value, 



by the steps of: (ci) ensuring that the received stabiliza- 
tion signal is within a preferred range of said known en- 
ergy value; (cii) calculating a value D„ that is proportion- 
al to a derivative of the spectrum at the stabilization peak 
energy, where n is the order of the derivative and is odd; 
(ciii) adjusting the gain of the tool on the basis of the 
value of D n . 
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